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Abstract 

Diseases affect coral species fitness and contribute significantly to the deterioration of coral reefs. The increase in frequency 
and severity of disease outbreaks has made evaluating and determining coral resistance a priority. Phylogenetic patterns in 
Immunity and disease can provide Important insight to how corals may respond to current and future environmental and/or 
biologically induced diseases. The purpose of this study was to determine if Immunity, number of diseases and disease 
prevalence show a phylogenetic signal among Caribbean corals. We characterized the constitutive levels of six distinct 
Innate immune traits in 14 Caribbean coral species and tested for the presence of a phylogenetic signal on each trait. 
Results Indicate that constitutive levels of some Individual immune related processes (i.e. melanin concentration, peroxidase 
and inhibition of bacterial growth), as well as their combination show a phylogenetic signal. Additionally, both the number 
of diseases affecting each species and disease prevalence (as measures of disease burden) show a significant phylogenetic 
signal. The phylogenetic signal of immune related processes, combined with estimates of species divergence times, 
Indicates that among the studied species, those belonging to older lineages tend to resist/fight infections better than more 
recently diverged coral lineages. This result, combined with the Increasing stressful conditions on corals In the Caribbean, 
suggest that future reefs in the region will likely be dominated by older lineages while modern species may face local 
population declines and/or geographic extinction. 
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Introduction 

Immune defenses are critical for species success on ecological 
and evolutionary time scales [1-3]. As species diverge, new sets of 
genetic, biological and/or environmental conditions are encoun- 
tered making it necessary for emerging species to trade off costs 
and benefits within and between traits [4], including those related 
to immunity [5] . Immunity plays an important role in the success 
of a given species and, in theory, evolves as species diverge [1,6], 
likely conserving beneficial mechanisms from ancestral species [7] . 
Depending on selective pressures (e.g. resources and/or stressors), 
the immune system develops novel strategies and diversifies during 
speciation [6,8], hence favoring individuals that survive patho- 
genic infections and other stressful events [7,9]. In closely related 
species, the study of immunity in relation with phylogeny can 
provide insight into the selective forces at work and the way 
organisms and populations may respond to them [6,10,1 1]. 

The deterioration of coral reef ecosystems has been associated, 
among other factors, with changes in environmental conditions 
(e.g. increased water temperature and ocean acidification) and a 
significant increase in number of coral diseases and epizootics [12- 
16]. However, the evolutionary importance of immune traits in 



corals has not yet been evaluated. In other organisms such as fleas 
[17], termites [18] birds [19,20], and vertebrates in general [21], 
immune traits are related to phylogeny. In corals, a recent 
taxonomic reorganization [22] facilitates assessment of trait 
variation and their relationship to life-history. 

Scleractinia, the Cnidarian Order grouping all reef-building 
corals, is divided into two divergent groups termed Robust and 
Complex corals [23-26], each composed of several non-mono- 
phyletic families [27] with dilferent evolutionary histories between 
the Atlantic and the Indo-Pacific regions [28] . Modern Caribbean 
scleractinians are grouped in at least six families, representing both 
Robust and Complex corals [27]. Recently, the Caribbean has 
become a disease 'hot spot' due to the high number of coral 
diseases, disease outbreaks and their widespread geographical 
range [16,29,30]. 

The innate immune system in corals is comprised of conserved 
components similar to those of other invertebrate [31,32] and 
vertebrate species [33,34] including the three general phases in the 
response to infection: recognition, signaling and effector responses 
[35]. While recognition receptors and several signaling pathways 
(e.g. toll and complement pathways) are activated upon pathogen 
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recognition, many immune components, such as some effector 
mechanisms, show constitutive activity (i.e. non-pathogen induced 
or basal levels). Some of the better studied immune mechanisms in 
corals include the melanin synthesis cascade (e.g. prophenolox- 
idase and phenoloxidase) [31,35,36], antimicrobial compounds 
[31,36-38], and antioxidants (e.g. superoxide dismutase, peroxi- 
dase and catalase) [31,32,35,36]. Combined, these immune 
components provide corals with the capability to control the 
presence and combat proliferation of pathogens [31], repair tissue 
[39], and reduce levels of reactive oxygen species generated during 
infections and associat(;d stress [40,41]. Constitutive levels of 
prophenoloxidase, melanin [35] and antimicrobial activity [42] 
have been linked to disease resistance [32,35,43-45] suggesting an 
active investment in components of the immune system [10]. 

Investment in immunity has been related to different life history 
traits [1]. Since all biological traits, including those involved in 
immunity, tend to vary within and across populations and/or 
species during speciation, different evolutionary pressures (e.g. new 
pathogens or changes in climate and/ or local environmental 
factors) can potentially result in distinct, species-specific, immune 
defenses. In corals, analyses of ecological strategies suggests related 
groups are affected similarly by environmental stressors and 
diseases [16,46-48], but litde is known about th(- role innate 
immunity plays in this pattern. Phylogenetic signals of biological 
traits are expected to be common across different groups of 
organisms [1 1,49], but it is unknown what types of traits or \vhat 
traits themselves will show a signal. Detecting a phylogt'nc-tic signal 
in coral immunity can provide insight into the evolution of the 
coral' immune system, and help explain the current pattern of 
disease resistance and its imphcations for the future of coral species 
and coral reefs, more so in fight of global climate change and 
increased disease pressure [1 1,50,51]. 

In this study, we tested for the presence of a phylogenetic signal 
in immune traits and epizootics in Caribbean corals. We 
characterized the constituti\c k'\'cls of six distinct immune traits 
in 14 of the most common and widely distributed coral species in 
the wider Caribbean. This represents ~20% of the total diversity 
of scleractinian corals in the region [52]. We also compiled 
published disease data and determined the levels of immunocom- 
petence of the studied sp(x:i(;s into two different mc'trics: Number 
of diseases (including both tissue loss diseases and growth 
anomalies) affecting each species, and mean prevalence, or the 
proportion of infected individuals in the population of a given 
species. These data sets were incorporated with the phylogeny of 
the host species (based on the 28S rDNA region) into three 
phylogenetic signal estimators (Bloomberg's K, Moran's I and 
Abouheif test). 

The null hypothesis was that in scleractinian corals there is no 
relationship between phylogeny and levels of various enzymes 
involved in constitutive immunity and no correlation with disease 
parameters. Our results show co-variation between constitutive 
immunity and phylogeny in Caribbean corals, with species in old(;r 
fineages (lineage defined as extant species and their ancestors) 
grouping together with lower number of disease and disease 
prevalence, than the modern lineages. Combined with other life 
history traits, these older lineages may be better suited to survive 
current and emerging diseases in the wider Caribbean. 

Materials and Methods 

Sample collection 

For this study, 14 of the most common Caribbean scleractinian 
species comprisinglO genera and 6 families [as defined by: 22, 27, 
28, 53], were collected from several reefs (Media Luna - 17° 



56.096 N; 67° 02.911 W, Turromote - 17° 56.097 N; 67° 01.130 
W, Conserva - 17° 57.831 N; 67° 02.940 W, Pinnacles - 17° 
55.963 N; 67° 00.714 W, Corral - 17° 56.986 N; 67° 00.504 W 
and Isla Cueva - 17° 57.599 N; 67° 04.827 W) off La Parguera, 
southwest coast of Puerto Rico (Table 1). To prevent seasonal or 
environmental effects, all samples were collected during the 
summer (northern hemisphere), the second week of August 
2012. This collection represents approximately ~20% of the total 
number of scleractinian corsd species in the region [54]. Many of 
these species, i.e. Montastraea cavernosa, Orhicella spp. ( = 
Monlaslraea annularis complex), Diploria, Pseudodiploria, and 
Poriles spp. are common and widely distributed through the 
region. Other groups, such as acroporids (A. palmata, A. 
cervicornis and A. prolifera) and pocilloporids (e.g. Madracis 
spp.) were not collected due to strict limits on sampling and 
manipulation. 

Small fragments from a total of 140 apparently healthy (i.e. with 
no signs of disease or bleaching) colonies (10 per species) were 
sampled. AU samples were collected under the specification of 
research collection permits to the Department of Marine Science 
University of Puerto Rico - Mayaguez (UPRM), issued by the 
Department of Natural Resources of Puerto Rico. A fragment of 
approximately 5 cm^ was carefully removed from the top of each 
massive/crustose colony with a hammer and a chisel. Small 
branches were broken from branching colonies. For the "free- 
living" Siderastrea radians, rolling stones larger than 4 cm in 
diameter were collected in shallow sea-grass beds next to the reefs. 
All samples were stored in indi\ idually labeled sterile Whirl-pack 
bags (Fisher Scientific, Waltham, MA), transported in seawater to 
the laboratory and flash-frozen in liquid nitrogen. Frozen samples 
were stored at — 80°C, shipped to the University of Texas at 
Arlington (UTA) in dry ice and kept at — 80°C until further 
analyses. 

DNA extractions and PCR amplifications 

The NCBI database has a significant number of sequences from 
most of the corals used in this project. The 28S rDNA region has 
been sequenced (as of September 2013) for 1 1 of the 14 species in 
this study (Table 1), and phylogenetic reconstructions showing 
similar topologies to other molecular markers, and divergence time 
estimation using fossils [Caryophyllia spp., Flabellum spp. and 
DendrophyMdae), were available. Sequences for three species 
{Meandrina jacksoni GenBank KJ946355, Pseudodiploria strigosa 
KJ946354 and Siderastrea radians KJ946356) were generated in 
this project after extraction of DNA using a modified protocol 
from Lajeunesse et al. [55]. A small fragment (~3 mg) of skeleton 
and tissue was mixed with (of glass beads (~200 |il, 1 mm, 
Ceroglass, Columbia, TN) and 600 ^1 of a cell lysis solution (0.2 M 
Tris, 2 mM EDTA, 0.7% SDS, pH 7.6) and shaken on a BioSpec 
(BartlesvUle, OK) beadbeater for 100 seconds. Proteinase K (3 (il 
— 20 mg/ml) was added and incubated at 65°C for 1 hour. The 
incubation was followed by protein precipitation with ammonium 
acetate (250 jtl —9 M) and freezing at — 20°C. The frozen extract 
was centrffuged (10,000 G for 15 minutes) and the supernatant 
removed, mixed with 600 (jl of isopropanol (100%) and centri- 
fuged (10,000 G for 5 minutes). The DNA pellet was washed with 
70% ethanol, air dried, and resuspended in 75 (tl of distilled water 
and stored at — 20°C. 

The 28S rDNA region was amplified using the 28SROM.IFw 
(5 '-GGCGACCCGCTGAATTCAAGCATAT-3 ') and 
28SDES.VRV 5'-GGTCTTTCGCCCCTATACTC-3') primers 
[56]. Reactions were performed using Perfect Taq Plus DNA 
Polymerase (5-Prime, Gaithersburg, MD) following the manufac- 
turer recommended reaction composition on 2 (tl of 1 :40 dilutions 
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Table 1. List of scleractinian coral species used to measure constitutive immunity and its variation across taxonomic levels. 



Family 


Genus 


Species 


NCBI accession 


Poritidae (III) 


Pontes 


Pontes astreoides 


EU262830 


Poritidae (III) 


Pontes 


Pontes pontes 


EU262878 


Siderastreidae (IX) 


Siderastrea 


Siderastrea radians^ 


KJ946356 


Meandrinidae (XII) 


Dendrogyra 


Dendrogyra cytindrus 


EU262819 


Meandrinidae (XII) 


Meandrina 


Meandrina jacl<soni* 


KJ 946355 


Meandrinidae (XII) 


Meandrina 


Meandrina meandrites 


EU262815 


Montastraeidae (XVI) 


Montastraea 


Montastraea cavernosa 


EU262810 


Merulinidae (XVII) 


Orbicetia 


Orbicella annularis 


HQ203479 


Merulinidae (XVII) 


Orbicella 


Orbicella franksi 


EU262849 


Merulinidae (XVII) 


Orbicelta 


Orbicella faveolata 


EU262781 


Mussidae (XXI) 


Diploria 


Diploria labyrinthiformis 


EU262772 


Mussidae (XXI) 


Pseudodiploria 


Pseudodiploria strigosa* 


KJ946354 


Mussidae (XXI) 


Mussa 


Mussa angulosa 


EU262869 


Mussidae (XXI) 


Mycetophytlia 


Mycetophyllia aliciae 


EU262809 


Underlined groups within the family and genus categories represent the sequence used for that group in the phylogenetic 
correspond to the sequence of the 28S rDNA region used for each species. (*= Species sequenced in this project). 


signal assessments. Accession numbers 



doi:l 0.1 371 /journal.pone.Ol 04787.t001 



of the extracted DNA (final reaction volume 25 |J.l). Amplifications 
consisted of 35 cycles of 95°C, 52°C and 72°C steps, each for 30 
seconds. Amplified products were cleaned with ExoSap (Affyme- 
trix, Santa Clara, CA) and sequenced with the forward primer 
using Big Dye 3.1 terminator mix (Applied Biosystems/Life 
Technologies, Grand Island, NY) on an ABI Hitachi 3730XL 
genetic analyzer at UTA' Genomics Core Facility. DNA sequence 
chromatograms were reviewed and edited using Geneious Pro 5.0 
[57]. The resulting sequences were combined with those obtained 
from the NCBI data base (657 to 685 bp) and alignments were 
performed on ClustalW using a gap-opening penalty of 1 5 and a 
extension penalty of 6 [58]. Phylogenies were constructed on 
MrBayes [59], using a general time-reversible model with gamma 
distributed rate heterogeneity (GTR+G) as substitution model, a 
chain length of 1,100,000 and 100,000 burn-in (phylogenies can 
be found in Text SI). 

In order to determine the approximate age of the studied 
lineages, that is the extant species and their ancestors, we 
calculated the divergence times on the 28S rDNA based phylogeny 
(Text S2). Divergence times were determined with BEAST 1.7.5 
[60], using a relaxed-clock uncorrelated lognormal allowing for 
nucleotide substitutions rates to vary between lineages [26]. The 
tree prior used a Yule process and the model of substitution was set 
to gamma distributed rate heterogeneity as suggested by 
jModeltest 2.1.3 [61], with invariant sites (GTR+G+I). Node ages 
and posterior probabilities were estimated on a run with 
10,000,000 generations and saving the topologies and parameters 
every 1,000 generations. In order to estimate the early divergence 
of the studied groups, calibrations were done using Dendrophylli- 
dae (~127Mya; including sequences of Tuhastrea coccinea, 
Cladopsamia gracilis, Leptosammia pruvoti, Endopachys grati, 
Enollopsamia rostrata and Balanophyllia spp.), Caryophyllia spp. 
(~160Mya) and Flabellum spp. (~77.5 Mya) as suggested by 
Stolarski et al. [26]. Phylogenetic reconstruction within BEAST 
was performed using Mean heights and node heights, a prior 
probability of 0.1 and 1,000,000 burn-in. 



Protein extractions and immune assays 

Protein extractions and enzymatic assays followed protocols 
previously used to study coral immunity [31,32,35,36,43]. Coral 
tissue was airbrushed from the skeleton using a Paasche single 
action artistic airbrush (Paasche Airbrush Company, Chicago, IL) 
with minimal amounts (5-6 ml) of Tris buffer (100 mM Tris, 
pH 7.8+0.5 mM dithiothreitol). To break open cells and extract 
proteins, tissue slurries were homogenized using a tissue homog- 
enizer (Powergen 125, Fisher Scientific, Waltham, MA) for 1 
minute on ice. One ml of the tissue slurry was added to pre- 
weighted 1.5 ml microfuge tubes for melanin concentration 
estimates. AH homogenates were centrifuged at 90 xG for 10 
minutes and the supernatant was recovered. Protein concentra- 
tions were estimated using the RED550 protein assay (G 
Biosciences, Saint Louis, MO) and standardized to a standard 
curve of bovine serum albumin. 

We performed assays for six immune traits: prophenoloxidase, 
melanin concentration, superoxide dismutase, peroxidase, catalase 
activity and inhibition of bacterial growth. Prophenoloxidase was 
tested on 20 |J,1 of the extract, mixed with 20 |j.l of sodium 
phosphate buffer (50 mM, pH 7.0) and 25 ^l of Trypsin (0.1 mg/ 
ml). The reaction was initiated by adding 30 jxl of dopamine 
(10 mM) as a substrate. Change in absorbance was measured 
every 30 seconds at 490 nm for 15 minutes and activity calculated 
during the linear range of the curve (1-5 minutes). 

Melanin concentration was assessed on the melanin-reserved 
portion of initial tissue slurry after freeze-drying (VirTis BTK, SP 
Scientific, Warminster, PA) for 24 hours. The resulting dried tissue 
was weighed and the melanin extracted with 400 |il NaOH 
(10 M). Extraction was done at room temperature for 48 hours at 
which time the tissue particles were centrifuged (90 xG) for 10 
minutes. 60 |J.l of the supernatant were used to determine the 
absorbance at 410 nm. Resulting values were standardized to a 
dose-response curve of commercial melanin (Sigma- Aldrich, Saint 
Louis, MO). 

Superoxide dismutase activity was determined with the SOD 
Determination kit (#19160, Sigma- Aldrich, Saint Louis, MO) 
following manufacturer's instructions. Absorbance at 450 nm was 
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measure in wells containing coral protein extracts and superoxide 
dismutase controls and compared to untreated samples. The 
inhibition was normalized by mg protein and presented as 
superoxide dismutase activity units per mg protein. 

Peroxidase activity was assessed on 10 |J,1 of extract with 40 |xl 
phosphate buffer (0.01 mM, pH 6.0) and 25 ^1 Guaiacol (25 mM). 
Activity was monitored for 15 minutes, recording the absorbance 
at 470 nm every 30 seconds. Peroxidase is presented as change in 
absorbance per mg protein per minute. 

Catalase was measured as the change in hydrogen peroxide 
concentration after mixing 5 |J,1 of the protein extract with 45 |il of 
sodium phosphate buffer (50 mM, pH 7.0) and 75 ill of 25 mM 
H2O2. Samples were loaded on UV transparent plates (Grainer 
Bio-one, Monroe, NC) and read at 240 nm every 30 seconds for 
15 minutes. Catalase activity was estimated as change in hydrogen 
peroxide concentration per mg of protein during the first minutes 
of the reaction. 

The percent inhibition of bacterial growth was assessed against 
Vibrio alginolyticus (Strain provided by K. Ritchie, Mote Marine 
Laboratory, GenBank # X744690). This particular bacterial 
strain was isolated from Orhicella faveolala, and has been 
implicated in Yellow Band Disease [62,63], and thus offers a 
good estimate of antibiotic activity in a broad range of corals. 
Bacteria were grown in salt amended (2.5% NaCl) Luria Broth 
(EMD Chemicals, Gibbson, NJ) for 24 hours prior to use in the 
assay. The resulting bacteria culture was diluted to a final optical 
density of 0.2 at 600 nm. 140 |al of the culture suspension was 
added to each well along with 60 |j.l of the coral extract. To detect 
possible effects of the media and the airbrushing buffer, controls 
with 60 |il of Tris Buffer (100 mM Tris, pH 7.8+0.5 mM 
dithiothreitol) were included on each experimental plate. Plates 
were incubated in the spectrophotometer for 6 hours at 29°C, 
determining the absorbance at 600 nm every 10 minutes. The 
change in absorbance during the logarithmic phase were used to 
determine the growth rate of each samples and the proportion of 
this rate to that of the bacterial control provided the percent 
inhibition for each samples. 

These immune assays are ideal for comparative immunity 
studies because they are not species-specific [10]. All assays were 
conducted on a Synergy 2 Microplate Reader (Biotek Instruments, 
Winooski, VT) and standardized to mg of protein when 
applicable. Constitutive levels of all six immune traits were 
compared between older and modern coral lineages (groups 
described in the results section). Comparisons were assessed with a 
t-test assuming unequal variances across traits. 

Data analysis and phylogenetic signal estimations 

In order to test the hypothesis that closely related taxa have 
similar activity of constitutive immune components, the data was 
partitioned in three groups corresponding to each of the clades of 
interest: family, genus and species (Table SI). Results from the 
immune assays were averaged for each taxon on each level. 
Additionally, to obtain an integrated measure of immunity the first 
component scores from a principal component analysis (PCI) with 
all individual immune measures was used as an additional 
category. Principal components reduce dimensions and convert 
multiple variables into composite indicators [64], improving the 
analysis of immune capacity in relation to other biological/ 
environmental parameters [65]. The analyses were done in JMP 
10.0.0 (SAS Institute, Gary, NC). 

In order to detect a phylogenetic signal, the immune data sets 
were compared against family, genera and species phylogenetic 
reconstructions (Text SI). The genus and family phylogenies were 
buUt with subsets of selected sequences from each group 



(underlined categories in Table 1). For example, Pseudodiplora 
strigosa 28S sequence was used as the representative for the 
Pseudodiplora genus. 

The possibility of phylogenetic signals in coral constitutive 
immune levels was tested with Bloomberg's K [49] and Moran's I 
measures as described by Gitdeman and Kot [66] and by 
Abouheif test [67]. K assumes the data follows a Brownian 
motion model (BM) and compares the observed phylogenetic 
signal with that of the trait under the BM model. Higher K values 
for a particular trait represent a stronger phylogenetic signal, and 
zero values indicate no effect of phylogeny [17,49]. Moran's I (7) 
on the other hand, is a model-independent measure of autocor- 
relation in which the relation between the variation in the trait and 
the phylogenetic distance is established. In this method, the data is 
divided into the phylogenetic component and the trait component 
and correlarograms are built to determine the effect of ranks and 
distances [66]. Lastly, the Abouheif [67] test (A), modifies Moran's 
I to successfully detect phylogenetic signal of different traits on 
phylogenies with both low and high number branches. All 
phylogenetic signal tests were performed on R using geiger, 
carper, picante, adephylo and phylobase packages. 

A Spearman's rank order correlation index was used to 
determine correlations between each immune measure and 
estimates of number of diseases and prevalence (proportion of 
infected individuals in the population of a given species at a given 
time) for each tsixonomic level (i.e. species, genus and family). 
Disease parameters (number of diseases and prevalence) were 
compiled from literature reporting epizootic events in the 
Caribbean from 1997 to 2005. Only manuscripts with species 
level resolution, and infected or diseased colonies specified as a 
subset of the total population surveyed were used so that number 
of diseases and prevalence estimates could be normalized across 
reports. Prevalence was estimated from the reports by combining 
data from all surveyed diseases (e.g. Black Band, White Band, 
White Plague, Yellow Band, Dark Spots, and growth anomalies) 
for each species. A total of 12 reports with appropriate information 
were found from different locations in the Caribbean, including 
Bermuda, Florida, Bahamas, Puerto Rico. Saint Croix, Bonaire, 
Yucatan-Mexico, Colombia and Venezuela (Table 2) [29,30,68- 
77]. The geographic coverage and extent of the disease surveys in 
these reports provided a very robust data set for both number of 
diseases and disease prevalence, thus truly representing the disease 
dynamics across the region. Phylogenetic signal in the number of 
diseases and prevalence was estimated as described above. 

Results 

Phylogenetic reconstruction and time of divergence 

Sequences of the 28S rDNA region for three corals were 

generated and in each case these newly sequenced species 
clustered according to current taxonomy; Meandrina jachsoni 
with its sister sp(;cies M. meandriles, Pseudodiploria slrigosa with 
Diploria labyrinthifomiis among the robust corals and Siderastrea 
radians with Porites spp. in the complex corals (Figure 1). The 
phylogenetic hypothesis obtained with the 28S rDNA region, 
using all 14 species, resolved the same groups delineated with 
other molecular markers and in other geographic locations 
[27,53,78]. 

Divergence times of the included lineages (i.e. extant species and 
ancestors) derived from the molecular clock show that Siderastrea 
radians (~247 Mya) and Porites spp. (~220 Mya) diverged early 
during the Triassic and appear as ancestral lineages among the 
Complex corals. Among the robust corals, the family Mean- 
drinidae also showed early divergence (~201 Mya), followed by a 
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Table 2. List of references, reports and reviews of coral disease parameters from various locations in the Caribbean used to obtain 
data for the 14 species used in the phylogenetic signal estimates. 



Reference 


Location 




Survey year(s) 


Data used 


Sutherland et at., 2004 


Caribbean Wide 


Review 


Number of diseases 


Weil, 2004 


Caribbean Wide 


Review of 1999-2002 


Number of diseases and prevalence 


Kaczmarski et oi, 2005 


Saint Croix 




2001 


Prevalence 


Voss and Richardson, 2006 


Lee Stocking Island - Bahamas 


2002, 2003 


Prevalence 


Santavi et oi, 2001 


Florida 




1997, 1998 


Number of diseases 


Ward et ai. 2006 


Yucatan Peninsula - Mexico 


2004 


Prevalence 


Garzon-Ferreira ef ai., 2001 


Colombia 




1998 


Number of diseases 


Navas-Camacho ef a\., 2010 


Colombia 




1998-2004 


Number of diseases 


Gil-Agudelo ef ai, 2010 


Colombia 




1998-2005 


Number of diseases 


Garcia ef ai. 2003 


Los Roques 


Venezuela 


1999 


Number of diseases and prevalence 


Croquer ef ai, 2003 


Los Roques 


Venezuela 


2000 


Number of diseases and prevalence 


Croquer ef ai, 2005 


Los Roques 


Venezuela 


2000, 2001 


Prevalence 



Data used included the number of diseases affecting each species and average disease prevalence (proportion of infected individuals on a population of a given 
species) for all diseases affecting each species. 
doi:1 0.1 371 /journal.pone.01 04787.t002 




COMPLEX 
CORALS 




Slderastreldae (IX) 



r 



strigosa 



D. labyrinthiformi^ 



1 



astreoides 



S. radians 



pontes 



Poritidae(lll) 





M. angulosa 
M. aliciae ^ 



D. cylindrus 
M. meandrites 

M. jacksoni 




Merulinidae 



O. annularis 

O. faveolata 
O. franksi 




ROBUST 
CORALS 



M. cavernosa 




eldae \ 
(XVI) 




Figure 1. Un-rooted Bayesian phylogenetic reconstruction (darl< line) using the 28S rDNA (657-685 bp) region from 14 species of 
Caribbean corals used in this study. The dotted lines represent the family groupings as recently proposed by Budd et al. [22]. All nodes in this 
tree have support (posterior probability) values higher than 0.80. 
doi:1 0.1 371/journal.pone.01 04787.g001 
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recent radiation into the current species (~ 30-42 Mya; Figure 2). 
The other Complex lineages evolved within the last ~100 to 
150 Mya (Figure 2). Faviinae (Z). lahyrinthifomis and P. strigosa) 
and Mussinae {Mycetophyllia aliciae and Mussa angulosa) 
diverged ~138 Mya, with current species diverging as recentiy 
as 40 Mya. The genera Orhicella and Montastraea originated 
~105 Mya, with Orhicella radiating between ~17 and ~7 Mya 
(Figures 1 and 2). This allowed the sampled lineages to be grouped 
into older (S. radians, P. porites, P. astreoides, D. cylindrus, M. 
meandrites and M . jacksoni) and modern (M. aliciae, M. angulosa, 
D. labyrinthiformis, P. strigosa, M. cavernosa, O. faveolata, 0. 
annularis and O. franksi) lineages. 

Phylogenetic signal in coral immunity and disease 

Results from the three phylogenetic signal measures were 
similar across taxonomic levels (i.e. species, genus and family). 
These analyses revealed different patterns of phylogenetic signal 
for each independent immune measure (Table 3). Significant 
correlations between phylogeny and immune traits were found in 
melanin concentration at the species level (Moran's I - /jp = 0.105, 
p = 0.036), peroxidase activity at the genus (/genus = 0.451, 
p = 0.006; ^g,,nus = 0.524, p = 0.005) and species levels (Bloom- 
berg's K - K,p= 1.269, p = 0.001; /,p = 0.488, p = 0.002; 
0.583, p = 0.002) and percent inhibition of bacterial growth 
family level (Kf^;iy= 1.637, p = 0.010, /fan^uy = 0.082, 



at the 



p = 0.048). Variation in prophenoloxidase, superoxide dismutase 
and catalase activities did not show a significant phylogenetic 
signal at any level (Table 3). 



Principal Component Analysis revealed that Principal Compo- 
nent 1 (PCI) explained 33.5% of the variation. Principal 
Component 2 (PC2) 25.1% and Principal Component 3 (PC3) 
23.3%. PCI values were similar between closely related species 
showing a significant phylogenetic signal at the genus 
(^ge„us = 0.325, p = 0.037) and species (/,p = 0.251, p = 0.036; 
yl,p = 0.271, p = 0.049; Table 3) levels. PC2 and PC3 did not 
show a significant phylogenetic signal. 

Closely related species have similar number of diseases and 
disease prevalence (proportion of individuals infected in the 
population of a given species). The number of diseases showed 
significant phylogenetic signal at the genus (^genus = 0.358, 
p = 0.037) and species (/:,p = 0.747, p = 0.015, /,p = 0.292, 
p = 0.026; y4sp = 0.351, p = 0.019) levels. Variation in disease 
prevalence was significant at the family (/family = 0.106, 



= 0.012) and species (K^p 



.102, 



= 0.009; /,p = 0.370, 



0.004; ^,p = 0.544, p = 0.007) levels (Table 3). 



Correlations between disease parameters and immune 
measures 

Several of the individual immune traits were significandy 
correlated to number of diseases and prevalence. Peroxidase 
(p = 0.971, p = 0.001) at the family level, prophenoloxidase 
(p = 0.667, p = 0.035), superoxide dismutase (p = 0.636, 
p = 0.048) at the genus level and prophenoloxidase (p = 0.564, 
p = 0.036) at the species level showed significant correlations with 
number of diseases. Only peroxidase (p = 0.971, p = 0.001) at the 



O. faveolata 




200 1 50 

Divergence time (Mya) 



Number of diseases 



Figure 2. Phylogenetic reconstruction with the 28S rDNA region used to estimate the divergence time (dark values in the nodes) of 
the evolutionary lineages of 1 4 Caribbean coral species. Time of divergence was estimated with BEAST following the method by Stolarski et 
al. [26]. The molecular clock was calibrated with specimens from the family Dendrophyllidae at 127±3.5 iVlya, and the genera Caryophyllia at 
160±3.5 Mya and Flabellum at 77.5±3.5 IVlya, represented by gray triangles and values. Number of diseases affecting each species is shown in the 
colored squares. Disease data was compiled from the literature (Table 2) and corresponds to surveys between 1997 and 2005 throughout the 
Caribbean. 

doi:1 0.1 371/journal.pone.01 04787.g002 
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family level was correlated with disease prevalence (Table 4). The 
integrated measure of immunity (represented by PCI) was 
positively and significantly correlated with number of diseases 
(p = 0.619, p = 0.018) and prevalence (p = 0.597, p = 0.024; 
Figure 3) at the species level. Number of diseases increases as 
divergence time of the lineages decrease. 

Older lineages such as Pontes spp. (known to be affected by at 
least 5 diseases), Siderastrea radians (5 diseases) and the Mean- 
drinidae species (4 diseases) are affected by fewer diseases than 
modern groups such as Diploria (7 diseases), Pseudodiploria (6 
diseases), Monlastraea (8 diseases) and Orhicella (8 diseases; 



Figure 2) [29,30,68-75,79]. Species from the subfamily Mussinae 
(Mycetophyllia aliciae and Mussa angulosd) were the exception 
with an average divergence time of ~33 Mya and a low 
susceptibility (3 and 1 diseases respectively). The fact that 
Mycetophyllia and Mussa form a well-supported monophyletic 
group [22,27,78] that according to our molecular clock diverged 
~138 Mya (Figure 2), may explain this inconsistency. Pairwise 
(t-tests) comparisons showed all six immune traits levels are 
significandy different (p<0.050) between older and modern 
lineages (Figure 4). 
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Figure 3. Correlation between (A) number of diseases (number of diseases affecting each coral species) and (B) prevalence 
(proportion of individuals infected in the population of a given species) with constitutive levels of immunity as measured with the 
Principal Component 1 (PC1) obtained from Principal Component Analysis of six immune measures (prophenoloxidase, melanin 
concentration, superoxide dismutase, peroxidase, catalase and inhibition of bacterial growth). The correlations were determined with 
the Spearman's rank order index (p) for both disease susceptibility and prevalence. Green squares represent older, and blue circles modern lineages. 
doi:1 0.1 371/journal.pone.01 04787.g003 
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Discussion 

Immune defenses are key to the evolutionary success of a species 
[3], but does the immune system retain ancestral characteristics, 
develop new or retain both ancestral and new defense strategies 
during speciation? [7] Answering this question in scleractinians is 
important since it can offer insight into the evolution of immunity 
and relevance to disease resistance and the future of corals in light 
of novel pathogens and climate change [50] . Results of this study 
indicate that the coral immune system, at least in Caribbean 
corals, has been molded by evolutionary history, with past selective 
pressures (biotic and abiotic) likely improving the innate immune 
system to respond more efficiently to new diseases [80]. Extant 
species, belonging to groups of species that have survived through 
more of these pressures, i.e. older coral lineages (Porites spp. 
~220 Mya, Sidemstrea radium ~247 Mya and Meandrinidae 
spp. ~201 Mya), seem to be evolutionarUy better equipped to 
cope/survive current conditions than more recent or modern 
lineages (Orbicella spp. = Monlastraea annularis complex, 
-105 Mya). 

Constitutive immunity, disease patterns and phylogenetic 
analyses of the 1 4 coral species in this study show some interesting 
patterns. Species from lineages that diverged more than 200 Mya 
are affected by fewer diseases, show lower disease prevalence, and 
have higher levels of some constitutive immune defenses. Some of 
these species, for example, P. astreoides and S. radians have 
ecological traits consistent with a resistant coral. Both are able to 
maintain normal physiological functions under stress [81], are 
brooders and exhibit a weedy-like dispersal and recruitment 
strategies [82]. Together with their high constitutive immune 
levels, these traits make for resistant and successful species [83] 
increasing survivorship and probably, species fitness. 

Some of the modern lineages comprise another distinct group. 
Species in this group are affected by up to 8 different diseases, 
show high disease prevalence levels and have lower constitutive 
levels of some immune defenses. Orbicella faveolata, O. annularis, 
Diploria lahyrinthiformis, and Pseudodiploria strigosa are among 
the species most heavily affected by disease in the Caribbean 
[16,29,30], leading to significant population losses in recent years. 

This study did not include the Caribbean Acroporids (Acropora 
palmata and A. cervicornis). Acroporids in general are thought to 
have very low levels of immunity [35] and the Caribbean species 
were nearly extirpated over their geographic distribution by White 
Band Disease [84,85]. High rates of clonal reproduction (i.e. low 
genetic variability) can result in higher population susceptibility to 
environmental and/or biological stressors [86]. Their clonal 
population structure [87] and hfe history traits (e.g. broadcasters 
with low recruitment success, fast growing branching morphology 
and high fragmentation rates, etc.) have rendered Caribbean 
acroporid populations highly susceptible to disease and contrib- 
uted to their significant losses over their entire geographic 
distribution. 

Overall, there was a relationship between constitutive levels of 
immune measures and diseases data that partitioned between old 
and modern lineages. Since the immune system is integrated, 
assessing several measures of immunity is ideal to gain a better 
understanding of both the potential to prevent and the capacity to 
fight infections [88,89]. This study did not include induced 
immune responses or resistance to pathogens in experimental 
settings. However, it appears our multivariate immune measures 
are good indicators of disease susceptibility across these 14 species. 
Similar patterns to those detected here have been reported for 
Indo-Pacific corals, where specific markers of the melanin 
synthesis cascade were related to disease prevalence [35]. Our 
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Figure 4. Proportion of tKie mean (± standard error) values of constitutive levels of six immune traits (percent inhibition of 
bacterial growth, melanin concentration, prophenoloxidase, superoxide dismutase, peroxidase and catalase activity) for older and 
modern Caribbean coral lineages (A). Pairwise comparisons for each trait using t-tests assuming unequal variances revealed significant 
differences between older and modern lineages for all traits. Older lineages showed higher mean values of melanin and percent inhibition of bacterial 
growth and lower values in the other traits compared to the modern lineages. Percent inhibition of bacterial growth and peroxidase activity (mean ± 
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standard error of change in absorbance 490 nm per mg protein per minute) for 14 Caribbean coral species (B). Black lines represent the average 
immune assay values for each coral lineage. 
doi:1 0.1 371 /journal.pone.01 04787.g004 



study integrates the melanin synthesis cascade, several antioxidants 
and general antimicrobial activity. 

Among the specific immune traits, inhibition of bacterial growth 
and melanin concentration were higher in the older coral lineages, 
suggesting that these two mechanisms (among others not described 
here) may be important in conferring resistance. Prophenoloxidase 
and all the antioxidants were higher in the modern lineages. There 
is an inverse relationship between the melanin product and 
proteins in its biosynthetic cascade. Prophenoloxidase activity is 
likely to increase upon microbial stimulation, but corals with 
higher melanin concentrations may be able to prevent infection by 
deploying a melanin barrier before the need to engage into 
additional, more costly responses [90]. 

All antioxidants were more active in the modern lineages than 
in the older coral lineages. Antioxidants are an integral part of the 
stress response to both abiotic and biotic stressors. AH modern 
lineages in this study are known to be susceptible to both elevated 
temperatures and to disease [31,36,43]. Increasing temperatures 
and/or pathogen load may keep these corals under a constant 
state of stress, thus elevating the levels of antioxidants at a cost to 
other defense pathways [32] . In addition to continuous stress, the 
algal symbiont type in each coral species may also play a part in 
oxidative stress that is reflected in the holobiont [91]. 

Variation in constitutive immunity among corals can also be the 
result of environmental and spatial gradients on a given reef 
While our data are from one region and one time point, we fell 
they are a start to developing hypotheses to why older corals may 
be more resistant to diseases. In fact, we attempted to reduce the 
effects of the environment by collecting during the same season 
(within a week), and either sampling the same species from 
different reefs and different species from the same reef and from 
the same location of the colony. The exception was Sideraslrea 
radians, which was collected in sea grass beds outside the reefs. 
However, S. radians immune measures clustered with those from 
both Porites porites and P. astreoides, supporting our hypothesis 
and minimizing the influenci; of the cm ironment. 

Associated organisms, Symhiodinium spp. and bacterial symbi- 
onts or communities may exert an important role in the immune 
defense of the coral host. Our protein extraction is enriched in host 
proteins but may include a minimal cimount of Symbiodinium 
protein. The corals in this study associate with a variety of 
Symbiodinium species/strains (clades A, B, C, D). At least seven of 
the 14 coral species (e.g. Pseudodiploria strigosa, Dendrogyra 
cylindrus, Meandrina meandriles, M. jacksoni, Orbicella faveolata, 
O. franksi and 0. annularis) are known to associate with 
Symbiodinium Bl [92]. The immune measures of these corals 
reflect the phylogeny of the coral host more than the presence of a 
particular symbiont. 

Caribbean reefs are not likely to disappear, but rather change in 
community composition and structure [,50,93] . Coral species with 
advantageous life history traits wiU be more likely to dominate [94] 
with the immune capacity of individuals and/ or populations 
serving as a critical piece of the puzzle. While our data are limited 
in the number of species and to the Caribbean region, the 
evidence presented here indicates that phylogenetically close 
scleractinians have comparable levels of constitutive immunity. 

Our results support the ongoing observations that disease 
resistance varies considerably among coral taxa and provides a 



possible mechanistic basis for this variability. Caribbean species 
from older lineages seem better adapted to resist infections with 
high levels of melanin and antimicrobial activity and, as a result, 
may have experienced significantly lower mortality from disease 
outbreaks in the recent past [95-97]. The life history traits of many 
of these older lineages seem to make them well suited for survival 
under the current changing conditions. Correspondingly, older 
Kneages of mammals appear less extinction-prone than previously 
predicted [98]. Recently diverged, modem Caribbean coral 
Uneages on the other hand, have been severely affected by 
diseases and/or thermal stress with significant tissue loss and 
colony mortalities at local and regional geographic scales and 
some of these taxa could face local extinction in many areas [99] . 

Supporting Information 

Table SI Mean values (± standard error) of six immune 
traits (prophenoloxidase, melanin, superoxide dismu- 
tase, peroxidase, catalase and percent inhibition of 
bacterial growth), the integrated immune measure 
(PCI) and disease parameters (number of diseases and 
prevalence) used in the phylogenetic signal analyses on 
14 Caribbean corals at three taxonomic levels (species, 
genus and family). 
(XLSX) 

Text SI Phylogenetic reconstructions based on the 28S 
rDNA region of 14 Caribbean corals at three different 
taxonomic levels (species, genus and family). These trees 
were used as input in the phylogenetic signal analyses. 
(NEWICK) 

Text S2 Sequence alignment of the 28S rDNA region 
used to determine the divergence times of 14 Caribbean 
corals. Dendrophyllidae (~127 Mya; including sequences of 
Tuhaslrea coccinea, Cladopsamia gracilis, Leptosammia pruvoti, 
Endopachys grati, Enollopsamia rostrata and Balanophyllia spp.), 
Caryophyllia spp. (~160 Mya) and Flabellum spp. (~77.5 Mya) 
were used as calibration point in the molecular clock. 
(NEX) 

Acl(nowledgments 

The authors would like to acknowledge the Department of Marine 
Scienees, University of Puerto Rieo Mayagiiez (UPRM) for partial funding 
for boat and diving, logistical support and laboratory' space. The authors 
would also like to thank Derek Soto, Duane Sanabria and Luis Rodriguez 
(UPRM) for their help with collection and sample manipulations, Whitney 
Mann and the stafTof the Genetics Core Facility (UTA) support in the Held 
and laborator\' work respecti\'ely. Comments by Matthew Fujita (UTA) 
and Jean-Francois Floi (MPFDS) and two anonymous reviewers helped 
improve this manuscript. The analyses were run on a server provided by 
the UTA Office of Information Technology, Samples were collected under 
the specification of research collection permits to the Department of 
Marine Science UPRM. 

Author Contributions 

Conceived and designed the experiments: JHPC LDM EW. Performed the 
experiments: JHPC EW JBT, Analyzed the data: JHPC, Contributed 
reagents/materials/analysis tools: JHPC LDM EW, Cfontributcd to the 
writing of the manuscript: JHPC LDM EWJB-L. Field work: JHPC EW. 



PLOS ONE I www.plosone.org 



11 



August 2014 I Volume 9 | Issue 8 | e104787 



Phylogenetic Signal in Coral Immune Traits 



References 

1. Prcvitali MA, Ostfcld RS, Keesing F, Jolles AE, Hanselmann R, et al. (2012) 
Relationship brtwcrn pace of life and immune responses in wild rodents. Mol 
Ecol Res 121: 1483-1492. 

2. Schulenburg H, Kurtz J, Morct Y, Siva-Jothy MT (2009) Introduction. 
Ecological immunology. Philosophical Transactions of the Royal Society of 
London, Series B: Biological Sciences 364: 3-14. 

3. Satake H, Sekiguchi T (2012) Toll-Like receptors of deuterostome invertebrates. 
Front Lnmunology 3: 1—7. 

4. Norris K, E\ an.s MR (2000) Ecological immunology: life history trade-offs and 
immune defense in birds. Behav Ecol 11: 19—26. 

5. Murphy PM (1993) Molecular mimicry and the generation of host defense 
protein diversity. Cell 72: 823-826. 

6. Loker ES, Adema CM, Zhang S-M, Kepler TB (2004) Invertebrate immune 
systems - not homogeneous, not simple, not well understood. Immunol Rev 198: 
10-24. 

7. Magor BG, Magor KE (2001) Evolution of effectors and receptors of innate 
immunity. Dev Comp Immunology 25: 651—682. 

8. TeUaJL, Scheuerlein A, Ricklefs RE (2002) Is ceU-mediatcd immuniu i claied to 
the evolution of life-history strategies in birds? Proc R Soc B 269: 1059—1066. 

9. Mcdzhitov R, Janeway CA (1997) Innate immunity: the virtues of a nonclonal 
system of recognition. Cell 91: 295-298. 

10. Millet S, Bennett J, Lee KA, Hau M, Klasing KG (2007) Quantilying and 
comparing constitutive immunity across avian species. Dev Comp Immunology 
31: 188-201. 

1 1 . Gilbert GS, Webb CO (2007) Phylogenetic signal in plant pathogen-host range. 

Proe Natl Acad Scl USA 104: 4979-4983. 

12. HarvcU D, Markel S, Jfirdan-Dahlgren E, Raimundo L, Rosembrrg E. et al. 
(2007) Coral disease, emironmrntal dri\'ers and the balance between coral and 
microbial associates. Oceanography 20: 36—59. 

13. HarvcU D, Altizer S, Gattadori IM, Harrington L, Weil E (2009) Climate change 
and wildlife diseases: When does the host matter the most? Ecology 90: 912—920. 

14. Weil E, Smith GW, Gil-Agudelo D (2006) Status and progress in coral reef 
disease research Dis Aquat Org 69: 1—7. 

15. Croquer A, Weil E (2009) Changes in Caribbean coral disease prevalence after 
the 2005 bleaching event. Dis Aquat Org 87: 33-43. 

16. Weil E, Rogers CS (201 1) Coral reefs diseases in the Atlanlic-Carilibean. In: Z 
Dubinsky and N Stamblcr, editors. Coral Reefs: an ecosy stem in transition. Dor- 
Dordrecht: Springer, pp. 465—491. 

17. Krasnov BR, Poulin R, Mouillot D (2011) Scalc-depcndencc of phylogenetic 
signal in ecological traits of ectoparasites. Ecography 34: 114—122. 

18. Buhner MS (2004) DupHcation and Diversifying Selection Among Termite 
Antifungal Peptides. Mol Biol Evol 21: 2256-2264. 

19. Lee KA, Wikelski M, Robinson WD, Robinson TR, Klasing KC (2008) 
Constitutive immune defences correlate with life-history variables in tropical 
birds. Mol Ecol Res 77: 356-363. 

20. Tieleman Bl, Williams JB, Ricklefs RE, Klasing KC (2005) Constitutive innate 
immunity is a component of the pace-of-life syndrome in tropical birds. 
Proceedings of the Royal Society, Series B Biological Sciences 272: 1715-1720. 

21. Richards MH, Nelson JL (2000) The evolution of vertebrate antigen receptors: a 
phylogenetic approach. Mol Biol Evol 17: 146-155. 

22. Budd AF, Fukami H, Smith ND, Knowlton N (2012) Taxonomic classification of 
the reef coral family Mussidae (Cnidaria: Anthozoa: Scleractinia). Mol Ecol Res 
166: 465-529. 

23. Romano S, Cairns S (2000) Molecular phylogenetic hypotheses for the evolution 
of scleractinian corals. Bull Mar Sci 67: 1043-1068. 

24. Romano SL, Palumbi SR (1996) Evolution of scleractinian corals inferred from 
molecular systematics. Science 271: 640-642. 

25. Romano SL, PEilumbi SR (1997) Moleculsu' evolution of a portion of the 
miitochondrial 16S ribosomal gene region in scleractinian corals. J Mol Evol 45: 
397^11. 

26. Stolarski J, Kitahara MV, MiUer DJ, Cairns SD, Mazur M, et al. (2011) The 

ancient evolutionary origins of Scleractinia revealed by azooxanthcUatc corals. 
BMC Evol Biol 11: 316^ 

27. Fukami H, Chen CA, Budd At\ Collins A, Wallace CC, et al. (2008) 
Mitochondrial and nuclear genes suggest that stony corals are monophyletic but 
most families of stony corals are not (Order Scleractinia, Class Anthozoa, 
Phylum Cnidaria). PLoS ONE 3: e3222. 

28. Fukami H, Budd A, Paulay G, Sole-Cava A, Chen CA, et al. (2004) 
Conventional taxonomy obscures deep divergence between Pacific and Atlantic 
corals. Nature 427: 832-835. 

29. Weil E (2004) Coral reef diseases in the wider Caribbean. In: E Rosenberg and Y 
Loya, editors. Coral health and disease. Berlin: Springer, pp. 35-67. 

30. Sutherland K, Porter JW, Torres C (2004) Disease and immunity in Caribbean 
and Indo-Paeifie zooxanthellate corals. Mar Ecol Prog Ser 266: 273-302. 

31. Mydlarz LD, McGinty ES, Harvell CD (2010) What are the physiological and 
immunological responses of cored to climate warming and disease? J Exp Biol 
213: 934-944. 

32. Palmer CV, A4cGinty ES, Cummings DJ, Smith SM, Bartels E, et al. (2011) 
Patterns of coral ecological immunology: variation in the responses of Caribbean 
corals to elevated temperature and a pathogen elicitor. J Exp Biol 214: 4240- 
4249. 



33. Hemmrich (i. Miller DJ, Bosch TCG (2007) The evolution of immunity: a low- 
life perspective. Trends Immunol 28: 449-454. 

34. Ghosh J, Lun CM, Majeske AJ, Saechi S, Schrankel CS, et al. (2011) 
Invertebrate immune diversity. Dev Comp Immunology 35: 959—974. 

35. Palmer CV, Bythell JC, Willis BL (2010) Levels of immunity parameters 
underpin bleaching and disease susceptibility of reef corals. The FASEB Journal 
24: 1935-1946. 

36. Mydlarz LD, Couch CD, Weil E, Smith GW, Harvell CD (2009) Immune 
defenses of healthy, bleached and diseased Montastraea faveolata during a 
natural bleaching event. Dis Aquat Org 87: 67—78. 

37. Vidal-Dupiol J, Ladriere O, Meistertzheim AL, Eoure L, Adjeroud M, et al. 
(2011) Physiological responses of the scleractinian coral Focillopora damicornis 
to bacterial stress from Vibrio (•oralliilyticus . ] Exp Biol 214: 1533-1545. 

38. Vidal-Dupiol J, Ladriere O, Dcstoumicux-Carzon D, Sauticre PE, Meistertz- 
heim AL, et al. (2011) Innate immune responses of a scleractinian coral to 
vibriosis. J Biol Chem 286: 22688-22698. 

39. Gonzalez-Santoyo I, Cordoba-Aguilar A (2011) Phenoloxidase: a key compo- 
nent of the insect immune system. Entomol Exp Appl 142: 1—16. 

40. Wang D, Li F, Chi Y, Xiang J (2011) Potential relationship among three 
antioxidant enzymes in eliminating hydrogen peroxide in penaeid shrimp. Cell 
Stress and Chaperones 17: 423-433. 

41. Bogdan C, RoUinghoff M, Diefenbaeh A (2000) Reactive oxygen and reactive 
nitrogen intermediates in innate and specific immunity. Curr Opin Immunol 12: 
64-76. 

42. Gochfeld DJ, Aeby GS (2008) Antibacterial chemical defenses in Hawaiian 
corals provide possible protection from disease. Mar Ecol Prog Ser 362: 119- 
128. 

43. Mydlarz LD, Palmer CV (2011) The presence of multiple phenoloxidases in 
C^aribbean recf-building corals. Comparative Biochemistry and Physiology, 
Part A 159: 372-378. 

44. Hawley DM, Altizer SM (2010) Disease ecfilog\' meets ecological immunologv: 
understanding the links between organismal immunity and infection dynamics in 
natural populations. Funct Ecol 25: 48-60. 

45. Boughton RK, Joop G, Armitage SAO (2011) Outdoor immunology: 
methodological considerations for ecologists. Funct Ecol 25: 81—100. 

46. Diaz M, Madin J (2010) Macroecological relationships between coral species' 
traits and disease potential. Coral Reefs 30: 73-84. 

47. Darling ES, Alvarez-FiUp L, 01i\ rr TA, McCUanahan TR, C.:6ic IM, et al. (2012) 
Evaluating life-history strategics of reef corals from species traits. Ecol Lett 15: 
1378- 1386. 

48. Roder C, Aril C, Bayer T, Aranda M, Daniels C, et al. (2014) Bacterial profiling 
of White Plague Disease in a comparative coral species framework. The ISME 
journal 8: 31—39. 

49. Blomberg SP, Garland T, Ives AR (2003) Testing for phylogenetic signal in 
comparative data: behavioral traits are more labile. Evolution 57: 717-745. 

50. Hughes TP, Baird D, Bellwood D, Card M, ConnoUy SR, et al. (2003) Climate 
change, human impacts, and the resilience of coral reefs. Science 301: 929-933. 

51. Roder C, Arif C, Daniels C, Weil E, Voolstra GR (2014) Bacterial Profiling of 
White Plague Disease across corals and oceans indicates a conserved and distinct 
disease microbiome. Mol Ecol 23: 965-974. 

52. Cairns SD (1999) Species richness or recent scleractinian. AtoU Res Bull 459: 1- 
45. 

53. Huang D, Licuanan WY, Baird AH, Fukami H (2011) Cleaning up die 
"igmessidae" Moleculsu" phylogeny of scleractinian corals from Faviidae, 
Merulinidae, Pectiniidae and Trachyphylliidae. BMC Evol Biol 11: 37. 

54. Weil E (2005) Current status of the marine biodiversity of Puerto Rico. In: P. M. 
a. E Klein, editor editors. Caribbean Marine Biodiversity: The known and 
unknown. Lancaster, PA USA.: DESteeh Publications Inc. pp. 85-108. 

55. Lajeunesse TC, Loh W, van Woesik R, Hoegh-(iuldberg O, Schmidt GW, et al. 
(2003) Low symbiont diversify in southern (ireat Barrier Reef Corals, relative to 
those of the Caribbean. Limnol Oecanogr 48: 2046—2054. 

56. Barbeitos M (2007) Phylogenetics and morphological evolution of Scleractinian 
corals. Biology. Buffalo: State University ol New York at Buffalo, pp. 125. 

57. Drummond A, Ashton B, Cheung M, Heled J, Kearse M, et al. (2009) Geneious 
v4.8. Available: http://wwwgeneiouscom/. 

58. Thompson J, Higgins D, Gibson T (1994) CLUSTAL W: improving the 
sensitivity of progressive multiple sequence alignment through sequence 
weighting, position-specific gap penalties and weight matrix choice. Nucleic 
Acids Res 22: 4673-4680. 

59. Huelsenbeck J, Ronquist F (2001) MrBayes: a program for the Bayesian 
inference of phylogeny. Bioinformaties 17: 754—755. 

60. Drummond AJ, Rambaut A (2012) BEAST is a cross-platform program for 
Bayesian MGMC analysis of molecular sequences. Available: http:/ / 
beastbioedacuk, 

61. Darriba D, Taboada GL, DoaUo R, Posada D (2012) jModelTest 2: more 
models, new heuristics and parallel computing. Nature Methods 9: 772-772. 

62. Gervino JM, Hayes RL, Poison SW, Poison SC, Goreau TJ, et al. (2004) 
Relationship of Vibrio species infection and elevated temperature to yellow 
blotch/band disease in Caribbean corals. Appl Environ Microbiol 70: 6855— 
6864. 



PLOS ONE I www.plosone.org 



12 



August 2014 I Volume 9 | Issue 8 | e104787 



Phylogenetic Signal in Coral Immune Traits 



63. CcrvinoJM, Thompson F, (iomcz-Gil B, Lorcncc E, (jorcau L'J, ct al. (2008) 
The Vibrio core group induces yellow band disease in Caribbean and Indo- 
Pacific reef-building corals. J Appl Microbiol 105: 1658-1671. 

64. Ma Z, Zhao J, Wu L, Shi C, Zhang G (2010) The quantitative research of 
composite immune indicator for crustacean. Fish Shellfish Immunol 28: 187- 
192. 

65. Hoi-Lcitner M, Romero-PujEinte M, Hoi H, Pavlova A (2001) Food availability 
and immune capacity in serin (Serinus serinus) nestlings. Behav Ecol Sociobiol 

49: 333-339. 

66. (iittlrman JL, Kot JM (1990) Adaptation: statistics and a null model for 
estimating phylogenetic effects. Syst Biol 39: 227-241. 

67. Abouheif E (1999) A method for testing the assumption of phylogenetic 
independence in comparative data. Evol Ecol Res 1: 895-909. 

68. WardJR, Rypicn KL, Bruno JF, HarvcU GD, Jordan-Dahlgrcn E, ct al. (2006) 
Goral diversity and disease in Mexico. Dis Aquat Org 69: 23—31. 

69. Voss JD, Richardson EE (2006) Goral diseases near Lec Storking Island, 
Bahamas: patterns and potential drivers. Dis Aquat Org G9\ 33 10. 

70. Kaezmarsky E, Draud M, Williams EH (2005) Is there any relationship between 
proximity to sewage effluent and the prevalence of Coral Disease? Caribb J Sci 
41: 124-137. 

71. Garzon-Ferreira J, Gil-Agudelo D, Barrios L, Zea S (2001) Stony coral diseases 
observed in southwestern Caribbean reefs. Hydrobiologia 460: 65-69. 

72. Garcia A, Croqucr A, Pauls S (2003) Estado actual dc las enfermedades y otros 
signos de dcterioro coralino en sietc arrecifes del Parque Nacional Archipielago 
dc Eos Roqucs, Venezuela. Rev Biol Trop 51: 173-180. 

73. Santa\y DE, Mueller E, Peters E, MacEaughlin E, Porter JW, ct al. (2001) 
Quantitative assessment of coral disease in the Florida Keys: strategy and 
methodology Hydrobiologia 460: 39-52. 

74. Navas-Camacho R, Gil-.^^udelo DE, Rodriguez-Ramirez A, Reyes-Nivia MC, 
Garz6n-Ferreira J (2010) Coral diseases and bleaching on Colombian Caribbean 
coral reefs. Rev Biol Trop 58: 95-106. 

75. GU-Agudelo DE, Navas-Gamacho R, Rodriguez- Ramirez A, Reyes-Nivia MC, 
Bejarano S, et al. (2010) Enfermedades coralinas y su investigacion en los 
arrccifes eolombianos. Bol Investig Mar Gosteras 38: 1—36. 

76. Groquer A, Pauls SM, Zubillaga AE (2003) White plague disease outbreak in a 
coral reefs at Eos Roques National Park, Venezuela. Rev Biol Trop 51 39-45. 

77. Cr6quer A, Weil E, ZubiUaga A, Marquez S (2005) Impact of a white plague-II 
outbreak on a coral reef in the archipelago Eos Roques National Park, 
Venezuela. Caribb J Sci 41: 815-823. 

78. Fukami H, Budd AF, Paulay G, Sole-Cava A, Chen CA, et al. (2004) 
Conventional taxonomv obscures deep divergence between Pacific and Atlantic 
corals. Nature 427: 832-835. 

79. Croqucr A, Bone D (2003) Eas enfermedades en corales eseleractinios: Un nuevo 
problema en el arrecife de Cayo Sombrero, Parque Nacional Morrocoy. Rev 
Biol Trop 51: 167-172. 

80. Veizer J, Godderis Y, Francois LM (2000) Evidence for decoupling of 
atmospheric C02 and global climate during the Phanerozoic eon. Nature 
408: 698-701. 

81. Muthinga N, Szmant AM (1987) The effects of salinity stress on the rates of 
aerobic respiration and photosynthesis in the hermatypic coral Siderastrea 
siderea. Biol Bull (Woods Hole) 173: 539-551. 



82. Harrison P (2011) Sexual reproduction of scleractinian corals. In: Z. Diibinsky 
and N. Stambler, editors. Goral reefs: an ecosystem in transition. Springer, pp. 
59-85. 

83. Lewis JB (1997) Abundance, distribution and partial mortality of Siderastrea 
siderea on degrading coral reefs at Barbados, West Indies. Mar PoUut BuU 34: 
622-627. 

84. Aronson RB, Precht WT (2001) White-band disease and the changing face of 
Caribbean coral reels. ll\drol)ifilogia 160: 25 38. 

85. Gladfelter WB (1982) White-band disease in Acropora palmata: implications for 
the structure and growth of shallow reefs. Bull Mar Sci 32: 693—643. 

86. Reusch T, Ehlers A, Hammerli A, Worm B (2005) Ecosystem recovery after 
climatic extremes enhanced by genotypic diversity. Proc Nad Acad Sci USA 
102: 2826-2831. 

87. Baums IB, Miller MW, Hellberg ME (2006) Geographic variation in clonal 
structure in a reef-building Caribbean coral, Acropora palmata. Ecol Monogr 76: 

503-519. 

88. Palmer GV, IVaylor-Knowlrs N (2012) 'I'owards an integrated network of coral 
immune mechanisms. Proe R Soc B 279: 4106—4114. 

89. Ardia DR, Gantz JE, Schneider BC, Strebel S (2012) Costs of immunity in 
insects: an induced immune response increases metabohc rate and decreases 
antimicrobial activity. Funct Ecol 23: 732-739. 

90. Christensen BM, Li J, Chen C-C, Nappi AJ (2005) Melanization immune 
responses in mosquito vectors. Trends Parasitol 21: 192—199. 

91. McC^inty ES, Pieczonka J, Mydlarz ED (2012) Variations in reactive oxygen 
release and antioxidant activity' in multiple Symhiodinium types in response to 
elevated temperature. Microb Ecol 64: 1000-1007. 

92. Finney JC, Pettay DT, Sampayo EM, Warner ME, Oxenford HA, et al. (2010) 
The telative significance of host-habitat, depth, and geography on the ecology, 
endemism, and speciation of coral endosymbionts in the genus Symbiodinium. 
Microb Ecol: 1-14. 

93. Hughes TP, Tanner J (2000) Recruitment failure, life histories, and long-term 
decline of Caribbean corals. Ecology 81: 2250-2263. 

94. Baker AG, Glynn PW, Riegl B (2008) Climate change and coral reef bleaching: 
An ecological assessment of long-term impacts, recovery trends and future 
outlook. Estuarinc, Coastal and Shelf Science 80: 435—471. 

95. Weil E, Croquer A, Urreiztieta I (2009) Temporal variability and impact of coral 
diseases and bleaching in La Parguera, Puerto Rico from 2003—2007. 
Caribb J Sci 45: 221-246. 

96. McClanahan T, Weil E, Cortes J, Baird AH, Ateweberhan M (2009) 
Consequences of coral bleaching for sessile organisms. In: M. van Oppen and 
J Eough, editors. Goral Bleaching: Patterns, Processes, Causes and Conse- 
quences Ecological Studies. Springer -Verlag. pp. 121-138. 

97. Miller J, MuUer E, Rogers G, Waara R, Atkinson A, et al. (2009) Coral disease 
following massive bleaching in 2005 causes G0"() decline in coral cover on reefs 
in die US Virgin Islands. Coral Reefs 28: 925-937. 

98. Verde Arregoitia LD, Blomberg SP, Fisher DO (20 1 3) Phylogenetic correlates of 
extinction risk in mammals: species in older lineages are not at greater risk. 
Proc R Soc B 280. 

99. Carpenter KE, Abrar M, Aeby G, Aronson RB, Banks S, et al. (2008) One- 
Third of Reef-BuUding Corals Face Elevated Extinction Risk from Climate 
ChEinge and Local Impacts. Science 321: 560-563. 



PLOS ONE I www.plosone.org 



13 



August 2014 I Volume 9 | Issue 8 | e104787 



